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ABSTRACT 


The primary purpose of this thesis was to provide a 
method of checking the geometry and element connectivity 
input data for two finite element programs, ADINA and 
SAP IV. This preprocessor will accept the ADINA or SAP IV 
data deck, with minor modifications, and generate a graphi- 
cal display of the finite element model. The display is 
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A. GENERAL 

The continued developments and advancements of the 
finite element method this last decade have provided 
greater machine capabilities than ever before. Here at 
the Naval Postgraduate School (NPS), two of the more 
favorable finite element programs are SAP IV [reference 1] 
for linear analysis and ADINA [reference 2] for nonlinear 
analysis. However, with the large amounts of numerical 
input/output data and automatic mesh generation, it is 
impractical to check and reduce this data without a graphi- 
cal representation. Data checking is divided into two 
areaS: preprocessing and postprocessing. 

I BrepeBeessing 

Preprocessing is the checking of the input data 

deck. Errors in a finite element program occur basically 
in two areas. First, how close is the mathematical model 
(omndary Condirerens, loading conditions, material proper- 
ties, etc.) to the real problem? Second, are numerical 
errors present, or did misinterpretation of instructions 
occur in data deck preparation? Of the second type, most 
common errors are found in the geometry and element connec- 
tivity data. Preprocessing includes the forming of a 
graphical representation of the finite element model on 


which geometry and element errors are easily detected. When 
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node and element numbering options are available, the graph 
aids in the physical interpretation of the output data. 
Preprocessing is not a foolproof method of eliminating 
errors, but it does provide a tremendous advantage to the 
user. Preprocessors may be incorporated into the data check 
mode of the finite element program. However, to modify a 
large and complex program is dangerous. This may not be 
the best approach. A safer method is to develop a pre- 
processor which will read the finite element program deck 
separately with minimum modifications to that deck. 
2. Postprocessing 

Though not as important as preprocessing, post- 
processing is extremely helpful in output analysis. Probably 
the most common and useful type of postprocessing is the 
contour plot. Appendix C of reference 4 lists a program 
developed to produce contour plots of stress data from finite 
element models. Contour plots can easily be adapted to a 
Prssystem, Due sberequires plotting the contours on 2D 
surfaces, a bit more complicated. Two methods of post- 
processing of displacements are the plotting of a deformed 
model or placing scaled vectors at the nodes. Іп the cases 
where the displacements are small, multiplication by a 
magnification factor produces an exaggerated representation. 
Like preprocessing, postprocessing can be incorporated into 
the finite element program directly, or done separately. 
When done separately, the finite element program must still 


Demmoatitea Slaghtly to obtain a punched deck of the stresses 
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and displacements in the desired format. This thesis is 


primarily concerned with preprocessing. 


B. HISTORY OF DEVELOPMENT OF PSAPl 

A package [reference 4] containing digital computer 
programs for generating oblique orthographic projections 
and contour plots was produced by the National Aeronautics 
and Space Administration's (NASA) Langley Research Center 
Perey ana distributed by the National Technical Information 
Service in January, 1975. The programs are completely 
general. Both programs contain options for selecting 
a i ous plotting equipment including CALCOMP, VARIAN, and 
cathode ray tube (CRT) displays. With minor modifications, 
they can be adapted to most any system. 

1. SUBROUTINE PSAP Implementation 

Losh [reference 6], for his master's thesis in 

aeronautical engineering, implemented the preprocessor and 
postprocessor program, PSAP, at NPS in December, 1976. 
Modifying the LRC package [reference 4], Losh adapted 
SUBROUTINE PSAP to the NPS IBM 360/67 system using the 
CALCOMP model 765 plotter. PSAP serves aS a preprocessor 
for SAP IV models, and serves as a postprocessor for dis- 
placements of those models. Unfortunately, PSAP is severely 
limited in the type of elements it can plot. 

2. Motivation for SUBROUTINE PSAPLI 

With the introduction of the ADINA [reference 2] 


Program atelrPowtn January, 1977, and with expectation of 
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doing future analysis on ceramic turbine blades, it was 
desired to expand PSAP to include all ADINA elements and 
the 8-20 node brick elements in SAP IV. Like PSAP, PSAP1 
contains preprocessing and displacement postprocessing 
capabilities.  PSAPl is presented in this thesis and has 
the following improvements over PSAP: 

a.  Preprocessing for all ADINA elements. 

b. SAP IV 8 and 8-20 node elements. 

c. Expansion of SUBROUTINE ERROR. 

ini ro Lion OF curves using shape functions 
[reference 3] through three points on the edges of the 
8-20 node brick elements and the 4-8 node plane elements. 

e. Improvements in defining the plot origin. 

II EIGN GO an option (ISCALE - 0) to plot 
sections of a model without losing perspective. 

g. Several other minor modifications considered 


improvements. 


CEC. PRESENT CAPABILITY 

PSAP1 has the capability to plot all ADINA elements 
and all SAP IV elements except the pipe element. It will 
interpolate curves on the edges of 4-8 node plane elements 
and 8-20 node brick elements. Many options are given in 
Appendix A. Some of the more frequently used options are 
listed below: 

i Numbering ef grid points (NOTAT = 1). 

2. Numbering of the elements (NOTAT - 2). 


DEP odd plot (KDISP - 2). 
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4. wPosbpaeecessing of daàsplacements (NUDISP or NVDISP 
or NWDISP = 1) in two forms: plot of deformed structure 
(KDISP = 1) or displacements represented by vectors at 
the nodes (KDISP = 3; see reference 6). 

5. Symmetric representation about the XY (KSYMXY = 1), 
XZ (KSYMXZ = 1) or YZ(KSYMYZ = 1) planes. 

Е Optaen to pet sections of themmodel (partial"plot) 
to obtain a better view. Partial plots may be plotted to 
the scale of the complete model to avoid losing perspective 
(ISCALE = 0) or blown up to obtain a better view (ISCALE = 1). 


Multiple plots may be obtained using the same geometry and 


same displacement data (KODE = 1), same geometry and new 
displacement data (KODE = 2), or new geometry and new 
displacement data (KODE = 3; see figure l). 


СЕЛЕСИ Е рле plots (sections, partial plots, 
additional problems) present no problem. Plotting package 
user courtesy dictates that no more than 5 plots be placed 
ο. ο CABCOMPUDIOtter at any one time. Also, if the plots 
contain many elements (especially 8-20 node elements), it 
PomoOssible собол out of space in the plotting data sets. 


When this happens, you will receive 
ERCOR AHC210TrT STAE;, ABEND CODE IS: SYSTEM OB37 SYSPLOT. 


ου "πο οσα 1S split the run into two jobs. т? 
the job must be run on one job (i.e.,a large number of 
elements in the model or an assembly drawing where tne 


τοι leple plots is the same), then SYSPLOT space 
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may be increased [references 7 and 8] by adding the card 
//GO.SYSPLOT DD SPACE=(CYL, (needed space) ) ,SYSOUT=C 

Hist prior to card 

//GO.FT10F001 DD UNIT = SYSDA 


in Appendix A. It would be wise to seek advice from a 
consultant in Ingersoll 146 if additional plotting space 


is required. 


D. EASE OF MODIFICATION 

Both PSAP and PSAP1 are written to maintain as much 
generality as possible for ease of expansion and modifica- 
tion. Several FORTRAN statements, variables and subroutines 
are not used. They were left purposely unchanged. Although 
PSAP1 specifically reads ADINA and SAP IV data, it can 
easily be expanded to include any geometry and element data 
format. Simply study the read-in and storage methods (see 
Section II), and construct appropriate subroutines to read 


any particular format. 
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II. PROGRAM ORGANIZATION AND DESCRIPTION OF OPERATION 


A. PSAPL FLOW CHART 

Figure l is a condensed flow chart of PSAPl. Probably 
the most important information given on this chart is the 
sequence in which the data cards, NAMELIST OPTION and 
NAMELIST PICT are read. Remember, when generating a 
sequence of plots, once a parameter has been defined, it 
retains that value until it is reassigned. Note that when 
KODE = l or 2, the original values of NAMELIST OPTION and 
NAMELIST PICT are retained until they are changed. However, 
when KODE = 3 a new title card, NAMELIST OPTION and a set 
of problem data are read. All variables in NAMELIST OPTION 
and NAMELIST PICT are assigned their default values. A 
new problem begins in this case. It is important that the 
last NAMELIST PICT to be read must contain the value of 


KODE = Q. 


B. NAMELISTS AND EULER ANGLES 
T. NAMELIST OPTION 

Description and default values of NAMELIST OPTION 
are given in Appendix A. Basically NAMELIST OPTION varia- 
bles pertain to the given problem: the number of nodes, 
geometry format, displacement format, space between plots, 
and paper size. Default vaiues for NAMELIST OPTION are 
set, and NAMELIST OPTION is read at the beginning of the 


problem. Assigned values will remain until exit from the 
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πού Fow chart 
of program PSAP1 


Open plotting data sets. 
Define the origin. 






Read and plot title card. 

Assign defaults to &OPTION 
and &PICT. 

Read &OPTION. 

Read and print geometry 

and element data. 









KUDISP=KVDISP= 


Read and print displace- 
ment data. 






Read &PICT. 
Scale the plot. 

Rotate original coordinates. 
Plot the elements. 






Hot ticle at top of graph. 


KODE=3 


Position the pen. 
Call PSTOP and close data sets. 
осор 







Ην 





program (KODE = 0) or a new set of problem data is read 
(KODE = 3). 
2. "NAMBLIST PICT 

Like NAMELIST OPTION, NAMELIST PICT variable des- 
criptions and default values are given in ЕЕ A. 
Basically NAMELIST PICT variables pertain to a given plot. 
One NAMELIST OPTION may apply to several successive plots, 
but each NAMELIST PICT defines a unique plot. That plot 
may include the whole model, part of the model and any 
options defined in NAMELIST PICT. NAMELIST PICT also speci- 
fies the viewing plane through the Euler angles (figure 18, 
Appendix A). 

a. Oblique Orthographic Projections (Euler Angles) 

Dos zample or an Oblique orthographic projection 

of a finite element model is given in figure 13 in Appendix 
A. The model can be viewed in any selected orientation. 
Euler angle transformations are used to specify orientation 
of the model to be projected. As described in reference 4, 
this transformation resolves the coordinate system of the 
model to a principal viewing plane (i.e., сл 08, Y 2.) 
ug nuch Che шр а їз to be plotted. Prior to rotation, 
the model coordinate system (X,Y,Z) is coincident with the 
coordinate system containing the viewing planes (X Σο, 20): 
The viewing planes are fixed, and the model is rotated about 
its model coordinate system. The rotations (,8,¢) of the 
body about the model axes (X,Y,2) are shown in figure 18, 


Appendix A. The NAMELIST PICT variables KHORZ (horizontal 
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БЕБУ ϱΙΕΕιειοαισαὶ axis), PSI (U), THETA (80), and 
PHI (¢) specify the viewing plane and Euler angles. The 
order of the Euler angle rotations must be PSI, THETA 


and then PHI. Mathematical transformations are: 


X X 
O 

yf = (A (Ag [AQ (Y 

Z Z 
О 


Cos w zen O 0 


ΙΒ i = Sin y Cos Y 0 
0 0 l 
COS © 0 Sin 6 
[λα] = 0 1 0 
-Sin 8 0 Cos 6 
1 0 0 
[А6] = 0 Cos ó -5іп ф 


0 Sin ф Cos 6 


SUBROUTINE ROTAT calculates the transformation matrices for 
every NAMELIST PICT, except if ISCALE = 0. ISCALE = 0 


directs the scale of the plot to be the same as that of the 
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EM —————:———-——rfÜÉÓóÓóÓóÓÜÉ——— ————— 


ae AAA A 


E CU DIGEL Should a rotation occur with ISCALE = 0, 
the plot width could exceed the paper width. 
be Scaling 

The safest scaling method is automatic scaling 
(ISCALE = 1). The user may specify a scale (ISCALE = 2) 
and the plot origin (XORGN,YORGN), but one must be careful 
Mee to run the plotting pen off the graph paper.  ISCALE = 0 
1s a very useful option. The plot will use the same scale 
as the previous plot. It is useful in an assembly graph 
where examination of a mesh in sections without losing 
perspective is desired. Example 3, Section III, illustrates 
the option ISCALE = 0. When ISCALE = 1 in a NAMELIST PICT 
defining a partial plot, a "blow-up" of that section is 
obtained.  ISCALE cannot be zero in the first NAMELIST PICT. 

c. Partial Plots 

To develop a partial plot, three methods of 
segregating elements exist: first, by the X, Y, and 2 
cutting planes; second, by node numbers, and, third, by 
element numbers. If a model has an area where the elements 
are relatively small, a "blow-up" may be desired. Choose a 
numbering scheme or coordinates to define the section to 
be segregated using one of the methods above. Example 3 
Peroure 14, Section III) uses X, Y, and Z cutting planes to 
define the partial plots. Example 4 (figure 17, Section IIT) 


uses element numbers to section the plots. 
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C. NODAL POINT (GEOMETRY) INFORMATION READ-IN 

Nodal point data is read in by the GEOMn subroutines 
(GEOM1, GEOM2, and GEOM9, see figure 2). Since SAP IV 
and ADINA data decks are similar, SUBROUTINEs GEOMl and 
GEOM9 are also very similar. They are both constructed to 
read and generate data in exactly the same way as ADINA and 
SAP IV. All data not needed by PSAP1 is disregarded and 
the nodal point data is stored in array 222 (figure 3). 
After studying storage array 222 and GEOMl (or GEOM9), 
a user moderately familiar with FORTRAN programming could 
easily construct a user supplied subroutine (GEOM2) to 


read the nodal point data in any desired format. 


D. ELEMENT (CONNECTIVITY) INFORMATION READ-IN 

After reading and storing the nodal point data, the 
element data must be read. The GEOMn subroutine will read 
the element control card (NPAR, references 1 and 2). 
SUBROUTINE ELTYPE (figure 2) calls the proper element sub- 
routine to read the element data specified on the element 
control card. If several groups of elements are to be 
read, the process is repeated until all of the element 
groups have been read. Although the nodal point data is 
stored in array ZZZ, the element connectivity is read and 


stored on device 10 (disk). 


E. DISPLACEMENT DATA READ-IN FOR DISPLACEMENT POSTPROCESSING 
Displacement data may be read in (figure 1) by SUBROUTINE 


DATA9 (KDATA=9) or SUBROUTINES DATAI1 or DATAS (user supplied, 
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KDATA=1 or 5). When read, displacement data is stored 

in the last three columns of array 222 (figure 3).  PSAP1 

can postprocess displacements for both ADINA and SAP IV. 

The difficulty comes in obtaining a punched deck of cards. 
Reference 6 gives a description of how to obtain a deck of 
cards for SAP IV in a format acceptable to SUBROUTINE DATA9. 
ADINA has no such provision. However, when preprocessing, 

the displacement data will be omitted (NUDISP=NVDISP=NWDISP=0) , 


and this step will be by-passed. 


ШО PLOTTING LOGIC 

SUBROUTINE PLOTX (figure 4) is the main plotting routine. 
Since the nodal point data is stored in array 777 апа the 
connectivity is stored on device 10, it is a Simple matter 
to read the connectivity from device 10 (one element at a 
time), and connect the nodes as they are defined in references 
1 апа 2. For example, NEND = number of nodes defining the 
connectivity of a single element, NUMEL = the element number 
and NODE(NEND) is the array containing the connectivity. 
Device 10 contains this information successively for each 


element. The statement 


READ (10) NEND,NUMEL, (NODE (I) , Il=1,NEND) 


will read the element connectivity to be plotted. The 
8-20 node brick (ADINA and SAP IV) and the 4-8 node plane 
elements (ADINA) may have 3 points defining each edge. If 


the midpoint node is defined, then isoparameteric 
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START 


Read elements individually to be 
plotted from device 10. 
READ (10) NEND ,NUMEL, (NODE (I) ,I=1 ,NEND) 


Check if element is in bounds ae 
tobe plotted 


TES 






NEND=2 NEND=4 NEND=8 NEND=12 OR 20 





Plot 2 Plot SAP IV РТО ADINA 
поде 8 node 8 node plane 
elements elements Briek and 20 node 
elements brick elements 
and SAP IV 
20 node brick 
elements 


NO Check if all elements have been 
read from device 10. 


XES 


RETURN 


Figure 4. Flow chart for PSAPl plotting subroutine, PLOTX. 
Era. ΡΕ Contains the connectivity of the 
element being plotted. NEND = the number of 
nodes defining the connectivity of a single 
element.  NUMEL - the element number. 
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shape functions (reference 3) are used to interpolate along 
the three-node edges. These shape functions are identical 

to those used by ADINA and SAP IV, so the geometry represented 
graphically is identical to the problem solved in ADINA and 
SAP IV. SUBROUTINE CURVE does the interpolation with the 


following equations: 


К = М1 * X TA * X ο X 


Em == NTR Y IN on ИШО Y 


NS (Ss => 1.0) / 2.0 


N2 = ЕО (S - 1.0) 


ο cse 1.0) / 2.0 


i 
E 
O 
| A 
(n) 
A 
E 
o 


PSAPl uses the NPS plotting package [reference 5] 


subroutines. 
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CII FSAPl SAMPLE PROBLEMS 


The following examples have been chosen to illustrate 
Some OL the most useful options of PSAP1. Prior to 
attempting to use PSAP1, the user should have the problem 
defined and the cards punched in the format of references 
l and 2. Appendix A of this thesis gives a complete des- 
cription of deck preparation for PSAP1 here at NPS. This 
section should prove helpful in the understanding and 


interpretation of the options presented in Appendix A. 


A. ADINA EXAMPLES 

Ire cd concrete Beam (example l, figure 5) 

This example was chosen because it illustrates 

the importance of the exploded plot when more than one 
element group is used. It is taken from the ADINA manual 
[reference lj. Figure 6 is a listing of the data cards 
as they are prepared for ADINA. Figure 7 indicates how 
that deck would be modified for PSAPl. Note: load cards 
are removed, NAMELIST OPTION and NAMELIST PICT are added, 
and the title to be plotted on the graph is added in figure 
7. Otherwise, figures 6 and 7 are the same. Figure 8, 


part (a), illustrates an undistorted (KDISP = 0) PSAPI plot 


with the nodes numbered (NOTAT = 1). Figure 8, part (b), 
shows the same mesh in an exploded form (KDISP = 2) with 
the elements numbered (NOTAT = 2). Note how the truss 


elements are visible in the exploded plot. 
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BEAM DIMENSIONS 


; E 


IS 25 à 
Bons» ie ¿tota qero 


FINITE ELEMENT IDEALIZATION 


MATERIAL PROPERTIES: 
ge. = 3740 psi 


σι = 458 psi Figure 5. Example 1, ADINA 
2 i truss and 8 node plane 

Py steel ~ 44000 psi elements, Reinforced Concrete 

E = 6100 ksi Beam (Given on page 84, 
concrete reference 2). 

y = 02 

Еб > 30000+k si 

E, = 300 ksi 
steel 
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Example 1, ADINA input deck listing, page 1 of 2. 
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Figure 6. 
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Figure 7. 


PSAP1 input đeck listing, page l of 2. 


Example 1, 
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2. Flat Plate With Hole (example 2, figure 9) 

This is a well known problem with which one can 
calculate the stress concentration on a hole in a plate 
under axial tension. Figure 10 is a listing of the PSAP1 
deck set-up. The mesh is composed of ADINA variable 4-8 
node plane elements. Notice on figure ll, parts (a) and 
(b), how the interpolating shape functions round off the 
3-node edges. Part (a) has the nodes numbered (NOTAT = 1). 
Part (b) has the elements numbered (NOTAT = 2), and illus- 
trates the use of the symmetry option (KSYMXY=KSYMXZ=1). 
The symmetry option enables one to obtain a picture of the 
complete plate even though the model only consisted of a 


quarter plate with proper boundary conditions. 


LU SAP IV EXAMPLES 

1. SAP IV Truss Problem (example 3, figure 12) 

Figure 13 is a listing of the PSAP1 data deck. 

Figure 14 indicates how multiple partial plots can be used 
to obtain a better representation of the model. Part (a) 
of figure 14 is the complete model. Part (b) is the left 
half (XXMAX = 50'), and part (c) is the right half 
(XXMIN = 50', XXMAX = 1.0E20'). Notice also that for 
Parts (b) and (c), ISCALE = 0, which means succeeding plots 
have the same scale as the first. Had ISCALE equaled l 
(blow-up), then the width of the half view would have been 
the same as the complete model. Figure 14 size is limited 


by the NAMELIST PICT variable PLOTSZ. 
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Figure 9. Example 2, flat plate with a hole in 
tension.  ADINA 4-8 node plane element. 


Thickness = 1 inch 

Young's modulus = 30.0 X πῶ psi 
Poisson's ratio = .333 

Total elements = 18 

Total nodes = 44 


Distributed Load = 
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PSAPI input deck listing, page 1l of 3. 


Example 2, 


Figure 10. 
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Example 2, PSAP1 input deck listing, page 2 of 3. 
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Eure ll. Example 2, PSAP1 output graphs, page 1l of 2. 


— — A p R FS 
` 


12 4 Ad 
| 
| 
| 
37 38 139 
| 
| 
28 ^q | | 
31 | 
E 2 | 
: B x 
22 C | 
Ш) 11 ñ x 140 
19 2 j | 
va 2a i | 
У. 3 N 36 
4 5 N ο \ 
ч. Е T \ 
oa > | 3S | 
N35 126 \ 
ζ 17 \ à 
E 18 2? [38 αἱ 


(a) Actual mesh, nodes numbered. 
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Hx contznued, page 2 of 2. 


18 L? m 18 


| 





(b) Symmetric representation, elements numbered. 
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NAMELIST PICT (LEFT HAND SIDE) 
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PSAP1 input deck listing, page 2 of 2. 


Example 3, 


Figure 13. 


ANO part of input deck: 





Figure 14. Example 3, PSAPL output graphs. 






(c) Right hand 
side. 


(b) Left hand 
side. 


(a) Complete Ser ruc eurea: 





_2. Cylindrical Bar With Spherical Hole (example 4, 
Figure 15) 


This mesh could have several uses. Two might be 
to calculate stress concentrations if the bar is under 
axial load or to calculate loading if the void is under 
pressure (i.e., dispersed nuclear fuel pellet). Figure 16 
1s a listing of the PSAP1 data deck. Figure 17, part (a), 
is a representation of the complete model. Parts (b), 
(c), (d), and (e) of figure 17 are partial plots of the 
total structure using options of the undeformed structure 
(KDISP - 0) with node numbering (NOTAT - 1), and the 


exploded plot (KDISP - 2) with element numbering (NOTAT - 2). 
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(c) Elements 1-10, elements numbered, exploded plot. 
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(d) Elements 11-20, nodes numbered. 
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(e) Elements 11-20, elements numbered, exploded plot. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


With a little study of PSAP1 and figure 2, one realizes 
how simple it would be to expand PSAP1 to include virtually 
an unlimited number of geometry data input formats. Simply 
supply the GEOMn subroutines, element subroutines and 
modify SUBROUTINE ELTYPE if necessary. Another interesting 
project would be to incorporate PSAPl into ADINA and SAP IV 
so that a mesh plot could be obtained in the data check 
mode only. Reference 4 contains routines for plotting stress 
contours which could be incorporated. Hidden line logic 
is available. It would be difficult but could be incorporated. 
Hidden line plots are pretty, but they are not very practical 
when preprocessing because in preprocessing it is advantageous 


to see all of the nodes. 
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APPENDIX A 


PSAP1 USERS MANUAL 


PAST TES To 


It may be useful for the user to review NAMELISTs in any 
good FORTRAN manual. However, a short description of the 
NAMELIST input is given here. Only columns 2 through 80 
of each card will be read. The computer expects to find a 
special delimiter symbol in column 2 of the first card 
followed by the NAMELIST name. The delimiter symbol for 
the IBM 360-370 series is the ampersand & ; other processors 
use the dollar sign $. Following the first card comes the 
list of variables and their assigned values. Some important 
points about the variables are: 

A. Variables and their values need not appear in any 
special order. 

B. Predefined variables need not appear in the list 
if it is not necessary to change the value. 

C. A comma should follow each assigned value. 

D. It is recommended that each variable appear on a 
Separate card to facilitate change. 

The last card following the NAMELIST variables contains, 
beginning in column 2, the delimiter sign followed by the 
Word END. For example, &OPTION (read NAMELIST OPTION) 


contains 12 variables. All are initially assigned default 
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values. But, suppose three of those values (NNDEST=200, 
KGEOM=9, YSPACE=2.0) do not apply to a specific problem. 
The &OPTION would be read as follows: 
_& OPTION 
_NNDEST=300, 
KGEOM=1, 


_YSPACE=5.0, 
| &END 


Blank space in column 1. 
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"ΤΗ PBEOCEDURES 


The following sequence of cards is necessary to use 


PSAP 1. 


NOTES IBM CARDS 
COLUMNS 


123456789-------------- 


(1) // [Standard green job card] 
// EXEC FORTCLGP,REGION.GO=150K 


A DORE NS YS IN DD * 


(2) e ΜΑΤΝ PROGRAM 
DIMENSION 222 (NZ) ‚DISPD(5,3,NON) 
CALL PSAP1(ZZZ,NZ,DISPD,NON) 
STOR 


END 


(3) ш 
//LINK.USDD DD UNIT=3330,VOL=SER=DISK02, 
// DISP=SHR,DSN=S1153.PSAPLl 
VIRENKASYSIN DD 7 
INCLUDE USDD (LOADM) 


ENTRY MAIN 
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(4) 


= 


(6) 


(8) 


(10) 


ОШ) 


(12) 


/* 
ἵ πο Ετιο ρο) DD UNIT=SYSDA, 
77 TSPACES (CYL, (3,1)), 


// DCB= (RECFM=VS , BULKSIZE=3520) 


//GO.SYSIN DD * 


[PSAP1 title card as it appears on plot] 


&OPTION 
[&OPTION variables to be initially set or changed] 


& END 


[ADINA or SAP IV geometry data. Title 
card through element data — remove load cards] 


[Case identification card — Omit if IDCASE = 0] 


[Displacement data cards — Omit for preprocessing 
only] 


&P ICT 
[&PICT variables to be initially set or changed] 


&END 


[Additional data as defined by last value of 
KODE in NAMELIST PICT — Omit if last 
value of KODE = 0] 


/* 
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Notes: 


(1) 


(2) 


eo) 


Standard basic deck set up as described in chapter 
3 of reference 3. 


// {Standard green job card] 
// EXEC FORTCLGP , REGION.GO=150K 
A BORT.  SYSIN DD κ 


Main program. 
DIMENSION 222(NZ) ,DISPD(5,3,NON) 
CALL PSAP1(ZZZ,NZ,DISPD,NON) 
STOP 
END 
The main program has two functions: to allocate 
fast storage space, and to call PSAP1. NON must 
be greater than the number of nodes. NZ must be 
+3 
greater than 4*NON(7*NON if displacement data cards 


are ES лр. i.e., NUDISP, NVDISP, or NWDISP - 1). 


Using load module library. 
/* 
//LINK.USDD DD UNIT=3330,VOL=SER=DISKO2, 
// DISP=SHR,DSN=S1153.PSAP1 
ο ΕΙ ο κ DD * 
INCLUDE USDD (LOADM) 
ENTRY MAIN 
PSAP1 should be precompiled and stored in the machine 
in a load module. Load modules are a type of user 
library described in reference 7 and chapter 3 of 


reference 8. Since PSAP1 requires over one minute 


to compile, precompiling results in a large time 
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(4) 


(5) 


(6) 


Savings. These cards are subject to change, and the 
most current version is determined by the user main- 
taining this library. If PSAP1 is not on a load 
module, these control cards may be replaced by 
subroutine PSAP1 (which consists of a box and a 


felt Of cards). 


Allocation of storage. 
/* 
//GO.FTLOFOO1 DD UNIT<SYSDA, 
// SPACE= (CYL, (3,1)), 
// DCB= (RECFM=VS , BULKSIZE=3520) 
PSAP1 uses a slow storage device to store the element 


connectivity. These cards allocate 3 cylinders as 


described in chapter 3 of reference 8. 
Deck set up card. 
// CO- SYSIN DD * 
This is a standard deck set up card described in 
chapter 3 of reference 8. 


PSAP1 title card. 


ОЛЕШЕН card is the title that will appear on the 
graph. Make sure a user identification is on this 
Ano consists of 80 alphanumeric characters. 

The first 40 characters will form the first title 


line. The last 40 will form the second line. 
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(7) NAMELIST OPTION. 


(Note: start in second column) 
δΟΡΤΤΟΝ 
[&OPTION variables to be initially set or changed] 
&END 
VARIABLE-DEFAULT DESCRIPTION 
ANDEST=200 Must be equal to the number of nodes. 
NUDISP-0 0 — X direction displacements not input. 
1 — X direction displacements input. 
NVDISP-0 O — Y direction displacements not input. 


1 — Y direction displacements input. 
NWDISP-0 0— 2 direction displacements not input. 
1 — Z direction displacements input. 
(Note: unless displacement data is to be input, 
allow NUDISP, NVDISP, and NWDISP to default.) 
KGEOM- 9 Specifies the geometry input format. 
l- Subroutine GEOM1 reads in ADINA data 
deck geometry and connectivity. 
2- Subroutine GEOM2 may be supplied by 
user along with subroutines NSPLAN, NS3DEE and NSTRUS 
to read nodal data and connectivity in any format. 
9- Subroutine GEOM9 reads in SAP IV 
data deck geometry and connectivity. 
KDATA-9 Specifies the subroutine and corresponding 
method of input for displacement data. 
1- Subroutine DATAl, a user supplied 


subroutine. 
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2- Subroutine DATA2, a user supplied 
Subroutine. 

9- Subroutine DATA9, reads a punched 
output displacement deck from execution of SAP IV as 
presented in reference 6. (Note: unless displacement 


data is to be plotted, allow to default.) 


NVALUS -0 Not incorporated, allow to default. 
IRESEQ-1 Not incorporated, allow to default. 
KPLOT-1 Not incorporated, allow to default. 
YSPACE-2.0 Space between plots in the Y direction 


in inches when successive plots are plotted (i.e., 
KODE # 0). The graph title is plotted both on the 
top and at the bottom of each set of graphs controlled 
by a given NAMELIST OPTION. The space between the 


title and the plot is YSPACE/2.0. 


PSIZE-9.0 Paper size in the X direction, in 
inches. Used in scaling of the plots to insure this 
dimension is not exceeded. However, when manual 


scaling (ISCALE = 2; see NAMELIST PICT) this protection 
is not available, and it is possible to exceed the 


paper width. 


TDCASE-0 0- No identification card precedes 


the deck of displacement values. 


1- Identification card precedes the 


deck of displacement values. 
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(8) 


(9) 


SAP IV or ADINA data cards. 


Here insert the geometry (node coordinates) and 
element connectivity. This includes the title card 
through the element data cards; the load cards are 
removed. For a SAP IV data deck, NAMELIST OPTION 
variable KGEOM = 9. For an ADINA data deck, 


KGEOM= 1. 


From ADINA deck remove: 
1. Applied loads data. 


2. Frequency and mode shape calculations data. 


From SAP IV deck remove: 
l. Concentrated load mass data. 
2. Element load multipliers. 


3. Dynamic analysis cards. 


Otherwise, these cards are exactly the same as the 
deck prepared for SAP IV or ADINA.  PSAPl is not 

limited to SAP IV and ADINA. The user may specify 
any unique format (i.e., KGEOM = 2). In this case, 


SUBROUTINE GEOM2 must be supplied by the user. 


Case ID card and displacement data cards. 


PSAP1 is intended to be used essentially for prepro- 
cessing. If it is desired to use the postprocessing 
option, the user is referred toreference 6. Otherwise, 


omit the case ID card and displacement data. 
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ο) ana (11) NAMELIST PICT. 


(Note: start in second column) 


&PTCT 


[&PICT variables to be initially set or changed] 
&END 


VARIABLE-DEFAULT DESCRIPTION 
KHORZ-1 Integer designating the horizontal axes 


PERT 2 


PHT-0.0 


DIETA-0. 


[Du 0. 0 


NEWFR-1 


ISCALE-1 


of the viewing plane. 1 - XC. 2 = Yo 3 = Zo ° 
(See fiçure 18.) 

Integer designating the vertical axes 
of the viewing plane. 1 = X ° 2 = Yo- 3 = Zo: 
(See figure 18.) 

Angular rotation of the model about its 
X axis in degrees (performed 3rd, see figure 18). 
0 Angular rotation of the model about its 
Y axis in degrees (performed 2nd, see figure 18). 

Angular rotation of the model about its 
Z axis in degrees (performed lst, see figure 18). 

1- Frame change before plotting. 

2- No frame change before plotting. 
(Normally allow to default. A frame change resets 
the Y origin past the previous plot by YSPACE given 
in NAMELIST OPTION and resets the X origin at 0.0) 

0- No scale change. Use the same scale 


as the previous plot. Useful in an assembly graph 


merecia 15 desired to examine a mesh in sections 


ES 
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without losing perspective. See example 3 in 
section III of this thesis. ISCALE cannot be zero 
in the first NAMELIST PICT. 

l- Automatic scaling of plot and 
computation of proper origin location. 


2- User specified origin and scaling. 


PLOTSZ-10.0 Maximum dimension desired on com- 
pleted plot prior to rotation. After rotation it 
is possible for the maximum dimension in the pro- 
jected plane to exceed PLOTSZ. The maximum size 
of the projection is unlimited on the vertical axis 
but is limited by the paper width (PSIZE in NAMELIST 
ОРОМ on the horizontal axis. If rotation of the 
model causes the projection to go off the paper, 
it is rescaled prior to plotting.  PLOTSZ is used 


ine c a Mng only LL ISCALE = 1). 


XORGN-0.0 Που στους Of Ele plot origin. 


(Used only af [SCALE = 2.) 


YORGN-0.0 loca ron tot the plot origin. 


(Used only if ISCALE = 2.) 


BSCALE-1,0 Model size reduction factor. 
PSCALE is equal to the actual model size diviaed by 
the desired plot size. it is used only if ISCALE = 2. 
(Note that when ISCALE = 2 is used, no rescaling 
occurs if rotation causes the projection to exceed 


the paper width.) 





NOTAT-0 O- No numbering on plots. 
Ἱ- Numbering of grid points. 


2- Numbering of elements. 


XLHT-0.14 Height in inches of the integers 
specified by NOTAT. It should be a multiple of 
ПЕСО XEHT IS not a multiple of .07, it will be 
rounded to the nearest multiple.  XLHT has a 


maximum of .49 and a minimum of .07. 


KDISP=O 0- Plot of undeformed structure. 

l- Plot of deformed structure. 

2- Exploded plot. 

3- Displacement represented by 
vectors. 

KDISP = l or 3 represents a form 
of postprocessing and displacement data must be 
ЕЭ). TF postprocessing is desired, refer 


to reference 6. 


IDMAG-2 l= Direct magnification of displace- 
ment data by DMAGS. 
2- Scaling of displacement data 


to a maximum value of DMAGS. 


BNAGS-1.0 Magnification of displacements 
ο τσ - ἹοςἈ. Reduction factor ot elements 
ADS? = 2). 

KO Y MX Y = 0 l- Symmetry about X-Y plane. 
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Roy MxXZ—0 l- Symmetry about X-Z plane. 


KSYMYZ-0 1- Symmetry about Y-Z plane. 
A plate quadrant with KSYMXZ and 
KSYMXY equal to l would yield a complete plate. 


See example 2 in Section III in this thesis. 


E To develop a partial plot, three methods of 
segregating elements exist: First, by X, Y and Z 
cutting planes; second, by node numbers; and third, 
by element numbers. The next ten variables are used 


to separate elements into partial plots. 


XXMAX , YYMAX , ZZMAX-1.0E20 Local cutting planes 


XXMIN,YYMIN,ZZMIN- (-1.0E20) Parallel to the principal 
planes. 


NDMAX-9999999 Maximum gridpoint identification 


number to be included in the plot. 


NDMIN-O Minimum gridpoint identification 


number to be included in the plot. 


NELMAX-9999999 Maximum element number to be 


mclwided in ἔπε plot. 


NELMIN-0 Minimum element number to be 


included in the plot. 
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KODE-0 


speeirfses the contre: option after 
the plot is completed. 

115; exit from program. 

1- Read another NAMELIST PICT. 

2- Read a new set of displacement 
data. (Postprocessing only.) For KODE = 2, dis- 
placement data must be followed by another NAMELIST 
PICT. 

3- Read a complete new set of input 
data starting with a title card. 

ROFMKODE = 1772 Or SS садатЕтопа! 
sections of the deck must be repeated. The deck 
must end with a NAMELIST PICT having a value of 


KODE = 0. 


Note: A most important point to remember when 
generating a sequence of plots is that once a 
parameter has been defined, it retains that value 
until it is reasSigned. For example, if PLOTSZ is 
assigned a value of 8.0 for the first of a series 
of plots, and it is not redefined in any subsequent 
NAMELIST PICT; the value of PLOTSZ will be retained 
as originally specified. However, when KODE = 3 
and a new title, NAMELIST OPTION and problem data 
are read, all variables in NAMELIST OPTION and 


NAMELIST PICT are assigned their default values. The 
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problem starts over in this case. Refer to the 
Flow chart in figure 1 and study the path for 


different values of KODE. 


(12) Delimiter card. 


/* 


Delimiter card is defined in chapter 3 of reference 


8. 





Ill. SPECIAL FEATURES OF PSAPI 


A.  POSTPROCESSING 

Reference 6 contains information and examples on the 
use of the postprocessor. SAP IV has the capability to 
punch displacement data cards in an acceptable format for 
PSAPI (specifically subroutine DATA9). As of this writing, 
ADINA does not have this capability, and data cards would 


have to be punched manually by the user. 


B. PARTIAL DATA 

Reference 6, on page 30, establishes a procedure by 
ΕΠΕ 15 possible to input only a portion of the finite 
element model for a data check. This feature is valuable 
in a case where several different people are preparing 
different parts of a large data base for a problem and 
desire to check individually their inputs graphically for 
accuracy. All nodal coordinates for the entire model may 
ИШ р 05 ОП those that specifically define the portion 
of the finite element model to be plotted. In either case, 
all nodal coordinates that define the elements to be plotted 
Base be speciried. To use this feature the "element control 
cards" (described in references 1 and 2) must be modified. 
For example, if only the connectivity for elements 15 
through 50 of a problem are available, the changes below 


would he made. 
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1. All SAP IV Elements 
All SAP IV element control card changes would 
have similar changes because columns 6-10 contain the 
total number of group elements. Columns 66-70 are not 
used. Thus to plot only elements 15 through 50, make the 
following changes: 


a. Enter the upper bound (1.e., 50) in columns 


Si Piter the lower Bound το. 15) in columns 
66-70. 
2. ADINA Truss, 2D and 3D Elements 
ADINA element control cards for the truss, 2D 
continuum and the 3D continuum elements would be changed 
as follows: 


em Enter the upper bound (i.e.; 50) in columns 


b. Enter the lower bound (i.e., 15) in columns 
9-56. 
3. ADINA Beam Element 
ADINA element control cards for the beam element 
would be changed as follows: 


a. Enter the upper bound (i.e., 50) in columns 


enter the lower bounayl1.e., 15) in columns 


653-68. 


Z 7 





MAIN PROGRAM 


Q OO OO 


PSAP1 


APPENDIX 
SUBROUTINE PSAP1 LISTING 


B 


сә с O O O O O O O O O OO O OQ O Q O O Q OO O O O O OQ O O OC O OC OQ 
HN OPN GP О С Осе ОО М оо О чемо М O O: O —4 CN cn +f uy O ο 
O oO oO O O oO O O Q rm ρήνιο ο ο οι οσο, 
©О о O O O O O O O O O O O O OCO CO O OC OO O OCO O O OC OQ OC OcCCOOOSO 
QOOOQOOQOQ]QQOQOoOoOOQ]Qonogoogooooogoogoooooooooooooco 
OOOOOOOOOOGOOoOoOOoOoOoooooooooooooooooooo 
ΟοοοοοοοοοοσοοοοσοοοοοοοοοσοοοσοοσοοοοσοοοώοσσοσσσοοΟ 
οοοοοοοοοωώωοοοοώοοώοσοοοοωώωοοώωοοώοώοοοοσοος 


SUBROUTINE PSAP1 DOCUMENTATION 


DESCRIPTION OF INPUT CATA CARDS 


O 
< O < 
ox 
.. a 
یاس‎ 
IU) 
< 
OL ν-- 
Ou Lu 
и E 
Z W 
I 

= ° 
ш ш 
= Z 
e Im 
"m Q. لے‎ 
F> <( 

a uu 
EO 
а a 
<и 
x Le 
¿Cç 

> 
u. #07 
Qoo 

-jm 
NULU 
σσ 
uy ш 
OT 
OZ 
«i «t 
x Z 
<í LU 
qp» 
WOW 

co 9 
Q«t-ul 
ч 2 
A Оа 
шш لے کک‎ 
ak 
ә SiS 
Zm Zz 


TITLE CARD = 


DED BY THE PROGRAM. 


TORAGE IN BLANK 
ΕΕ 


VERIFY 5 
VALUES N 


N AND CONTROL 


NAMEL IST OPTION - CONTAINS VALUES TO 
COMMO 
THE FOLLOWING VALUES ARE INCLUDED--- 


= 
MO 
Drm 
5 С 
4 ° ° 9 
ш = =< т^ 
DU СО Ce Q 
© = = = 
< == = Rh 
> x ο O ο) 
ЫЈ ш ш ЫЈ 
со α a x 
e = 9 =4 وسم‎ 9 94 
g5 ZO = ес 
LUZ Ot ol QI! 
Uv» > o FP 
әш = > = 
< Oe oz ¿(= 
иш بللا‎ We LL py 
с Се а 4 ce 
O un ке) жч) 
O< Or” Or Ok 
= | < | 1 2 
ЫЈ жи Sl NU 
NI = = = 
= тш ZW ш 
= sO MO "05 
(O < <£ < 
Or < <= <. 
а. а. Ба ra 
- IN un {am 
O Om Om Cem 
< 2) O e С) 
сс С? = = = 
OW = (9 292,209 
LUZ uz uz 
uc Zin Ze Ze 
OO WO WO WO 
OO DUDO 
ασ x. d ato 
ul $ JO JO JO 
ο O.Z* Q. Z *%* Q. > # 
== Q n уе k 
= On = к= 
σας CNMCO<Q< OC <( OC< O 
E = e + 
Li 1 O< II O< II OC !! 
p< κο ο GC 


Xu οἱ-- = = > 
E LN JLX A X md X X d 
mXICOFPF2SOODSOCOOCOO2 


= «ug «uus ud «t 
o» Uude tu. u. LL. LL. 
uy MO WO WO mU O mU 
0.0 O O a 
it H I H H H Hl 
+ 3t Χ t 
E O ж ва Χα. ἃ 
2 Y) Y) VY) 
LL) =) ο وسم‎ 
O e e (=) 
< =) > = 
= = 25 = 


СЕОМ1 


SUBROUTINE 


BROUTINE AND CORRESPONDING METHOD OF INPUT 
"SUPPLIED 


SU 
RY 
ER 


QOO OOQOOQCQ OOQ(OO(OQOQOQOQOOOQOQOQOQO OQO Oc UOOQO OVO OU 


78 








QOOUQOOOOOOOOQOOOOOOOOOOOOOOQOOQ]QQQ]Q]QoOoOooOoooogqdoooooocdooc 
DOAN APD DOM DOA OANMO FIND ОС Оч УСО О С О чс Оо СО Ус TINO 
ی ب کج ب د‎ ST SSAA OLA LAUD A MA LA DO о оро 9 OD 9 O Mo Po foo PA ММГ со со CO cO CO CO cO 
ΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟ 
οοοοοοοοοοσοοοοώοοώωοοοοουοοσοοοοώωοσοοοοοοοοοοοοσοσσσ 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOCOOCOOOOOOOOOOooooocooo 
οοοοοοοοοοσώοοοοοοοσώωοοοοοοοοοοοοώοσοοοοοοοσοσσσσοσο 
οοσοοοοοοοοοοοοοσοοοοωώοοώοοοοσοςοοοοοοοοοοοοο 


po 
m ~ 
~ m 
a D ч ч 
< ra rr BP 
= < 2 < «x 
ez, " ο ο 
I o 
= u E = 
<= OO O c 
= шш ш ш 
<< ddm OZ 2 ° 
O U. C ш uj < 
to I о о tr 
> O e x < < 
PAS ш а.з а 
кч Z a a 
Ш 5 ка N mn > 
Zorel OF = — m 
шш = r су > 
WITT m с -- 
OA O Z< < < > 
¡OT Z = w WY < 
SO © LL 
Z= uJ) QUO с с uj 
© > “( < < < O 
Q ш ш ш ш 
αι CO c c с = 
20W Œ 9 O 
EE Cc © roca = 
<A mO Ὁ Бшш 2 
wo ш ш y < 
К Ох) б ew2g 
NQ Z Da <4 l 
О ш> = > Piük | 
لان سإ‎ «X Ud «t IL. «T | 
ш w OTO- O | 
Omaa Z = | 
ш. к< = ш Wow 


aa e Z>ZOZ UY 
Oa #D De mat OW 
-- ww Oder FOr*Z 
UNDO CFF DOL a 
> maoan 
шс ш — C OO ar <1 
оша. си IVN оса с I 
wa eDa DaD o 
"DINO RNZNOA NIN au) 
= Z July 2 YN о 
OXXOIMHZXNX wDD 
шсошашшо CIOS 
шш оше ош {Чч 1 шш 


bad 
< O < 
e = 
< < 
OOO 
YU ا‎ 


NVALUS 


LL 9 
οι O 
LL «t uw 
1ο. 050 
° з? uL 
e COW 
ш г С) = 
VW) Y= „4 >С 
= ш „u 
> το <t 
لے‎ ш O ° Or 
> OR ZN WO 
< 3 4 <. 2. 
Le za زلا س‎ = 
wk CO oO eU) 
O = < = 
ш ш Qa. O 
O O uo pay Y) [4.27 ° 9 
Hp m= 2 > > WO uu 
> Q ı De Ow ш 
=> ш æ Us e (y ZN 
=) i= mila ZU OAOA 
M Ez OOO OZ w>z> 
< D W = at Dot 
Qa ο “о0о PO Uu O 
( = oz Ὁ o cuc 
= as «oc« uv) aWowW 
μ΄. ΒΗ O. m =uj >= 
I с > x ος ΟΙ 
| ш ej NWO Cr LOAD 
O u>zZ rodó qa 4 
U) NO O= XU OSJO 
C. uy wz لے‎ Mx Y + adja. 
Z а Ou AN a 2X2* UDN0aAU 
+ > ۾‎ N FU ev) J (C) Jf 
**=*p»0a0>uwuwI23*20-30 ZON OF 
«1 zx UU) We οἱ Opg LL 
о «шошо чиа. ва гм ои Ju O 


O ο. & 5 bk "0 Cr- Uu 
NN JOD I! = uJ ANO m I 
N <£ < Á Wj ў чок? 
HF-—FUIO(UOOLÉJZA Pour x XE 
- ου لے‎ кеш We IC ed 
MR Dae de охо ош SF DACDOWOW 
< Ə<( LL. O O C)O <( CO —JF— <( <( _J <í Cu OA 


LL. Z LL. = LL. LL. LL. L. Ü. << LJ cou u. 
U Wo WAS Ow wo es Uj 
O OUmM NNT OV Сл @ a 
а. 11 l! 1 
+С ПОП ПП! + + + 
+ WwW 3t + is + 
N Б uy Lu UN 
ш OO ы ον) < 
Co. καὶ} «t =ч Q 
= AA а. VY) O 
>< >< 2 Q. =$ 
> 


OOOOOVOOOHOODOOOODODOOODOLDODOHOODDOOOOVODODOOOO 


79 


RMS, 
ST OPTION. 


2 
OUTINE GEOM2 WHICH IS PREPAREO BY THE USER TO 


H 
1 


I 


F 
p 
UTINE GECM1 WHICH READS ADINA GEOMETRY DATA 


Y 
0 
М 
O 


OS Ч абд 


К 

Ü 

R 
GEOP 
UBR 


E 

Ν 

Γ 

L 
E IF KGE 
CALL SUB 


DE 

D 
US 
US 


WOUUO QOUOUOO 





ωοοωοσοοοοοοοοοοοοοοοοοοοοοοοοοωώωοοοοώώοσοσσοσοσώσο 
мос О CONO) HF UO pr 000 Су ιο TIn Or OOO Чоу ОМ с О ч ССОО М O O О cum? 
ωσωωσσσοσοσσσοσοσοσσοσοσςΟ)κ κκ κ ee e mM κκ ος NAIQUN AIAN CIM 0 00 0 0 
OO OOOO OO MOO OO ete m m mE e اسم ہے اس‎ md us tm uci emi emi e ed ad em rm unl d rm d md md e t uL e rd e ml 
ωοοωοοοοοοοοοώοώοοοοοοοοοοοοοοοωοοοοοοοώώοσοοουοσσοσἒς 
QOQOOOOOOQOOoOOOOOOOoOgOoooogocoooooooooooooooooooooocococ 
QOOOQOOOQQOQQQOOQOqQQoQouog]QOogOogooouoouooOog]ooOOoo]uooOoOodoooocooodcoooodoc 
οοωώωοοοοοοώοοοοώοσοοοσοοώωοοοοοοοώὄοσοοώωοώοωώωσοοοσώοσος 


READ GECMETRY DATA. 


= 9 


USE IF KGEOM 


GEOMETRY DATA. 


CALL SUBROUTINE GECM9 WHICH READS SAP IV 


CASE IDENTIFICATION CARD. 


а. шұ 


cur 


e 
UJ «t 


И. σΌοσ-α 
eu سم‎ 
C) X < 
uJ e «to 
OW) Ud 
„Duo. 
Oa zZ 
N «Xu 
— > > 
bul 
O OU 


uv) - 
NZ jo. 
To IY) 
O <A port 
OF =o 


Q Ou 
LLC) 00 => 


О.с 
Qao < 
WM 
> Dr 
Þe C2 
NO O 
zo F 
COT a 


{== == р: 
= >С 


pe 
O Ze zZ 
œx U j O 
«tU» 
WE ov 

e C <t 
na Ou 
om |... 
TUZO 
که را سم‎ 


NG FORMS, 
ELIST OPTION. 


FOLLOWI 
IN NAM 


1 


OUTINE DATA1 WHICH IS PREPARED BY THE USER 


TA 
R 


UTINE DATA5 WHICH IS PREPARED BY THE USER 


Ξ 9 


USE IF KDATA 


NAMELIST PICT - CONTAINS VALUES NEEDED TO GENERATE PLOTS. 
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